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The indirect carrier-phonon vertex function is evaluated for cubic zirconia by employing the many- 
body pertnrbation theory on top of the frozen phonon model. The phonon-induced carrier velocity 
and spectral renormalization are represented at each individnal lattice vibronic mode. Our results 
show that for the indirect electron-phonon scattering, optical phonons make up the dominant contri¬ 
bution. It is observed that, despite the phonon-phonon interaction, there is a substantial polaronic 
character for all the activated lattice optical vibrations which is in accordance with the previous 
experimental and theoretical finding concerning the crucial role of optical phonons. Eventually our 
results are capable to discuss that, in the presence of optical phonons a transition between small to 
large polaron is probable. 

PACS numbers: 


I. INTRODUCTION 

In semiconductors, insulators, and semi-metals with a 
small electron density of states at the Fermi level, lat¬ 
tice vibronic states (phonon) are the dominant quasi¬ 
particles that are responsible for the electronic and ther¬ 
mal conductivity. Accordingly, studying the phonon- 
induced effects in materials has attracted immense atten¬ 
tion. These interactions are identified to determine the 
highest operation limits of any electronic device. Strik¬ 
ingly, high dielectric materials are new emerging can¬ 
didates for improving and modifying the electronic and 
thermal conductivity in the new generation of electronic 
devices [Il-Q. Among a wide variety of these materials, 
zirconia(.Zr 02 ) is of significant scientific and technolog¬ 
ical importance due to its superior mechanical stabili^ 
and outstanding electronic and thermal properties [3-0 ■ 
One of the encouraging applications of Zr02 is its po¬ 
tential ability to be replaced by Si 02 , which is a conven¬ 
tional gate oxide of field effect transistors. Zr 02 has a 
high dielectric constant (e ^ 25) with a large band gap(5- 
6eV). Compared with the conventional silicon-based gate 
oxides, Zr 02 prevents leakage current more effectively [3] . 
There are many theoretical and experimental studies in 
electronic structure and lattice vibrations of different 
phases of Zr02^ l3-[l^. However, in contrast to the 
electron-electron(e-e) coupling, there are only a few avail¬ 
able theoretical and experimental studies on the electron- 
phonon(e-ph) interaction in Zr02. For instance, an ab 
initio investigation of the phonon dispersion and the 
lattice contribution to the dielectric properties of zir¬ 
conia is reported in Ref. [l3| . The calculated Born ef¬ 
fective charge in their study indicates the crucial role 
of the phonon vibrations to the static dielectric con¬ 
stant €o, which is higher for the cubic phase in com¬ 
parison to the other zirconia phases. In Ref. 0, El, 
using density functional perturbation theory(DFPT), a 


Longitudinal-Transverse optical(LO-TO) mode splitting 
at F point which is qualitatively large, is reported. This 
effect is a result of the non-analytic term that requires 
the knowledge of the Born effective charge tensors and 
dielectric constant of Zr02. According to their study, 
when taking phonon-phonon (ph-ph) coupling into ac¬ 
count, the phonon softening is inevitable and the re¬ 
ported Gruneisen Parameter also reflects the requirement 
of anharmonic approximation for the precise calculation 
of electronic and thermal transport. In this regard, the 
variation of phonon distribution should change the elec¬ 
tronic distribution, effectively. All the aforementioned 
arguments declare relevant clues about strong e-ph inter¬ 
action in zirconia. Recently, the electronic correlation of 
Zr 02 in the framework of many-body perturbation the¬ 
ory in GW approximation [i 3-[13| is explored in Ref. . 
Their results address the interplay between the volume 
deformation potential(hydrostatic strain) and the elec¬ 
tronic degree of freedom. They discussed the fact that, 
the dielectric permittivity tensor of zirconia is dominated 
by the static contribution(eo) while, the electronic part 
(coo) is a factor 5-6 times smaller than cq. Such a big 
discrepancy results in the emergence of a considerable 
electric polarization, particularly due to the LO phonons. 
This phenomenon highlights the e-ph coupling, which is 
neglected in the study of Jiang et ai. In Ref.[2l|, the 
obtained results from an ab initio calculations for dielec¬ 
tric function showed a good accordance with the experi¬ 
mental results, since they considered the electron-optical 
phonon contribution to evaluate the dielectric function 
of Zr02. Despite all difficulties, there are numerous ex¬ 
periments on the electronic structure and phase transi¬ 
tion properties of Zr02 compounds. The importance of 
the phonon scattering and its influence on the thermal 
conductivity of the Zr 02 microstructure is discussed in 
Ref. [1^ , where the change in the surface strain was mea¬ 
sured. In another experimental study using ellipsome- 
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try, the obtained dielectric function of polycrystalline Yt- 
tria stabilized zirconia(YSZ) thin films illustrated conse¬ 
quential electron coupling with the optical phonons 
Nevertheless, a comprehensive study that focuses on e- 
ph scattering effects on Zr 02 electronic features, is still 
lacking and the present work may compensate it for 
future studies in this field. 

In this study, we measure the vibronic states of cubic 
Zr 02 through indirect phonon activation in the frame¬ 
work of the frozen phonon model. Our work has the ad¬ 
vantage of being compared and scaled with those applica¬ 
tions in which the e-ph interaction plays a consequential 
role in tuning the efficiency of electronic devices through 
managing quantities such as the figure of merit. This pa¬ 
per is structured as follows: In Sec. II we will explain the 
theoretical approach, i.e. a brief theory of all quantities 
which are discussed in this study, and details of our tech¬ 
nical computing. In Sec. Ill we will reveal and discuss 
our findings with the inclusion of the combined many- 
body and DFPT, henceforth denoted GqW^DFPT and 
emphasize the electronic and thermal features of zirconia 
focusing on its cubic lattice phase. Section IV summa¬ 
rizes the main points and findings. 


II. COMPUTATIONAL METHOD 


Within the state-of-the-art density functional pertur¬ 
bation theory approach [l^. the direct e-ph vertex func¬ 
tions using the Migdal approximationis: 


Knik,q) = 



( 1 ) 

where 'ipn,k and ipn^k+q are the electronic Bloch states, 
Uscf is the self-consistent Kohn-Sham potential created 
by all electrons and ions, are the eigen vec¬ 

tor components of the dynamical matrix correspond to 
the displacement of i-th ion(mass Mi) in the direc¬ 
tion j due to the phonon mode (A, g), and Uq^ij = 
exp(iq.Ri)uiij, where uuj denotes the displace¬ 
ment of i-th ion from unit cell 1 in the direction j and 
Ri{l = 1,2,...,V) denotes the unit cell. Since the di¬ 
rect calculations of the e-ph coupling matrix using the 
aforementioned method requires lots of time consuming 
computations, we apply many-body perturbation theory 
in the GW approximation [ij, 1251 - 1^ by employing the 
displacement snapshots which are deriyed from DFPT 
in the framework of the frozen phonon picture. Thus 
GoWo@DFPT is capable of measuring the strength of 
the indirect e-ph coupling and carrier scattering in cubic 
Zr02- Having earned the e-e correlation in each indi¬ 
vidual snapshot lead us to the qualitative prediction of 
the coupling between the electronic cloud and phonon de¬ 
grees of freedom. Phonon branches comprise nine modes 
arranged in three longitudinal and six transverse parts for 
acoustic and optical phonons. Treating each phonon as a 
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FIG. 1: The phonon induced electron-electron interaction( 
phonon emission(a) and phonon absorption(b)). 


perturbation of the self-consistent Kohn-Sham potential 
is crucial to our theoretical study. 

Fig [1] represents the indirect e-ph interaction schemat¬ 
ically, one electron of wave vector (k) emits a phonon of 
wave vector (q) and scatters into the state of (k — q), then 
a phonon of wave vector (q) immediately absorbs by an¬ 
other electron of wave vector (k) that scatters into the 
state {k'+ g)®. Refer to FigH] the associated electronic 
states and their related energies are: 

|z) = \k,k'-0) 

|a) = \k-q,k';l) 

I/) = \k - q,k' + q;0) 

Ei = Ek -f Ek' 

Fa — Ej^ — q E}^r fiUJq 

Rf — Rk — q F Rk'+q 

assuming fiojq Huid where wd is Debye frequency, it is 
deduced that within the energy shell of the order of hujo 
around Fermi energy, this process leads to an attractive 
interaction between electrons. The indirect e-ph coupling 
strength is given by: 


{f\Rindirect\‘^} — \^q 


huJn 


\Ek- Rk-q)^ -MuJ^q 

TlUJn 


+ 


{Ek'+q-Rk')‘^-k^ujl 


( 2 ) 


where, the scattering amplitude for electron transition 
from the state (k) into the state (k-q) with the cre¬ 
ation of a phonon (q) on the vacuum phonon state is 
Mq{l\aq^\0) = Mq, similarly the scattering amplitude 
from the state (fc') into the state {k' + g) with the ab¬ 
sorption of a phonon (q) is: M*{0\aq\l) = M* 

In many-body calculations, to investigate the elec¬ 
tronic properties of the systems, the quasi-particle self¬ 
energy has the leading contribution. In this scheme, the 
carrier lifetimes or equivalently damping rates can be ex¬ 
tracted from the imaginary part of the self-energy IsTI- 

H: 


r 


-1 


-2 


dr 


dr'ipi{r)Im'E{r, r';ei)'ipi{r') 


( 3 ) 


where ipi{r) and Ei are the electronic states and the re¬ 
lated energies, respectively. Accordin g to Hedin[^ along 
with the GW on-shell approximation!^ the self-energy 
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expands at the first order in the dynamically screened 
interaction: 


S(r,r' : Ei) = ^ J de'e : e'-Ei)M^(r, r' : e') 

(4) 

in which, G is the Green function of the electronic sys¬ 
tem, W is the dynamically screened Goulomb interaction 
and (5^ denotes a positive infinitesimal time. In this ap¬ 
proximation, by neglecting the off-diagonal terms in the 
self-energy equation one can determine the poles of the 
Green function that result in the lifetimes of the elec¬ 
tronic excitations: 

AEij(w) = J dr J dr'tlj*{r)[T,{r,r',uj) - VLDAir) 

6{r - r')]'ijji{r') (5) 

in this equation vlda denotes the local one electron po¬ 
tential, and the electron damping rate follows as: 

T~^ = -2ZiImAT.ii{£i), ( 6 ) 


The quasi-particle spectral residue or energy renormal¬ 
ization factor(Z-factor) is an important quantity to eval¬ 
uate the e-ph coupling strength in the interacting many- 
body systems and is defined as Zn(k) = \{'^nk\iJnk)\‘^, 
where lipnk) is the bare electron state in the n-th Bloch 
band and is that of the perturbed system. The 

Z-factor characterizes the e-ph interaction strength and 
assesses how much quasi-particle energies are close to the 
DFT values. Moreover, in the case of no correlations the 
Z-factor is equal to the unity and generally defined as: 


Z^ = {l 


dReAIlii{uj) 

duj 


O—Si } 


(7) 


The quantity — is a parameter which de¬ 

notes the strength of the e-ph coupling. Thus a lower 
Z-factor indicates a stronger e-ph coupling. 

Another important feature of our many-body calcula¬ 
tions is the velocity renormalization factor which corre¬ 
sponds with the real part of the self-energy and is related 
to the e-ph coupling strength through: 


1 


- 


-1 


Vr{k) - Vo{k) 
Vi{k) 


( 8 ) 


in this equation Vi(k) and vo{k) denote the electron ve¬ 
locity in the non-interacting system and in the presence 
of e-ph interaction, respectively. 

In following we will emphasize our computational de¬ 
tails to achieve the e-ph coupling which causes the renor¬ 
malization effects and eventually results in polaronic fea¬ 
tures in our system. To perform the frozen phonon cal¬ 
culations in the GoWo@DFPT framework, we used the 
QUANTUM ESPRESSO package [^. and to improve the 
quasi-particle properties, the dielectric matrix was calcu¬ 
lated using the general plasmon-pole approximation as 


implemented in the Yambo package (3^. Valence elec¬ 
tronic wavefunctions are expanded in a plane-wave ba¬ 
sis that is truncated at a cutoff energy of 50 Ry. The 
electronic and phonon calculations were performed on 
a uniform 6x6x6 momentum grid, and the ultra-soft 
Martins-Troullier norm-conserving pseudo-potential with 
Perdew-Wang LDA exchange and correlation were ap¬ 
plied for both Zirconium and Oxygen atoms. 


III. RESULTS AND DISCUSSIONS 
A. frozen phonon calculations 

FigEl shows the GqWq@DFPT results for zirconia 
electronic density of states in which the unperturbed 
DOS is located in the background of each panel. The 
acoustic branches (a), (b), and (c), declare no deviation 
from the unperturbed DOS. In the other words, due to 
the presence of acoustic phonons, ions displaced in-phase. 
Hence, acoustic phonons have no effect on the energy re¬ 
distribution of the electrons. To put it another way, the 
interactions between acoustic modes and electrons are 
qualitatively energy conservative and the correlation be¬ 
tween the electrons as well as the coherency of the sys¬ 
tem, being maintained. As a consequence, to random¬ 
ize the electron velocity due to the elastic scattering, is 
the only effect of electron-acoustic phonon coupling[36l]. 
Panels (d), (e), (f), (g), (h), and (i) of Fig l2] depict the 
density of states for the six perturbed lattice vibrations 
in the presence of optical phonons. The electron-optical 
phonon scattering is an inelastic interaction; therefore, 
the electrons relax by emitting a LO or TO phonon as 
is illustrated in FigH] These inelastic scatterings cause 
considerable redistribution of the carrier momenta and 
energies as is shown in Figj^l Zirconia is an ionic ma¬ 
terial; therefore, the Frohlich interaction is expected to 
play the main role in the scattering processes. This inter¬ 
action is proportional to the [ojLoi^^ ~ Eq 
it it stated before, the ratio of cq and Coo varies between 5 
and 6 orders of magnitude, that confirms the relativity of 
the scattering processes to the presence of LO phonons. 
Apparently in optical phonon branches, the density of 
state satellite character is destroyed drastically from the 
perfect background, except for the panels (e) and (f). In 
these two modes, the band-gap of Zr 02 is still clear, but 
smaller than the perfect DOS; therefore, Zr02 behaves 
like a narrow band-gap semiconductor. 

FigEl shows the Z-factor for the conduction states in 
cubic Zr02- In the acoustic modes no discrepancy is ob¬ 
served compared to the background as is observed in the 
DOS pattern; therefore, no appreciable momentum re¬ 
distribution occurs during the acoustic phonon induced 
e-e interaction. Two results are depicted from the op¬ 
tical modes: 1) the second TO and the first LO modes 
in panels (e) and (f) exhibit an increase in energy renor¬ 
malization, contrary to the other optical modes, and 2) 
these two branches maintain the symmetry of the crystal 
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FIG. 2: (color online)Cubic zirconia density of states ver¬ 
sus electron energy. Energy redistribution occurs due to the 
phonon induced e-e interaction. 


FIG. 3: (color online) Renormalization factor due to the 
phonon induced e-e coupling. Renormalization in optical 


field. The Z-factor is a measure of correlation between an 
individual electron with the rest of the electronic cloud. 
As it is argued in Ref.fl^, the typical Z-factor values of 
the transition metal oxides vary from 0.5 to 0.8 for d and 
s orbitals and it is well known that a higher Z-factor re¬ 
sults in a lower correlation between electrons. From this 
point of view, compared to the other optical modes, it 
can be elucidated from snapshots (e) and (f) that, the 
e-ph coupling is reduced, indicating a higher values for 
electronic mobility. The calculated Z-factor for (e) and 
(f) show that, the phonon-induced electron renormaliza¬ 
tion is insufficient for these electrons to have polaronic 
character. In other optical modes, based on the com¬ 
parison between binding energy and the zirconia band¬ 
width together, with a reduction in the Z-factor, result 
in the polaronic nature of electronic transport in cubic 
Zr02- The considerable reduction in the renormaliza¬ 
tion weight (Z-factor), confirms the obtained results in 
the Ref.j^ in which the dielectric matrix of Zr 02 is re¬ 
produced by taking the phonon self-energy along with the 
electron self-energy, into account and also improvement 
in the obtained screening length by taking e-ph coupling 
into account, is in raod agreement with the available ex¬ 
perimental values 21|. Inferred from our calculations, it 
is also predictable that, increasing the Z-factor value in 
panels (e) and (f), and then decreasing drastically in (g), 
(h), and (i), results in a small to large polaron transi¬ 
tion by the substantial role of the optical phonons as is 
clarified in Ref. . The second effect of those phonons 
which break the crystal symmetry and remove the de¬ 
generacy of cubic Zr 02 , is that electrons in our system 
to behave like electron gas with lacking the single quasi¬ 
particle character and simultaneously following the Fermi 
liquid model of quadratic pattern for the scattering rates 

Ref. [ 3 . 

The calculated velocity renormalization of electrons is 
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FIG. 4: (color online) Velocity renormalization versus electron 
energy. Polaronic and non-polaronic features for carrier trans¬ 
port is deduced from the velocity renormalization curves as is 
described in the text. 


reported in FigU) In acoustic branches we observe no 
important changes in the velocity renormalization factor 
which means that the activation of the acoustic modes 
during in-phase ionic displacements should not affect the 
momentum redistributions of electrons. Furthermore, 
the out-phase ionic displacements in (e) and (f), result 
in the velocity renormalization promotion, confirming 
a non-polaronic nature in the first LO and second TO 
snapshots, consistent with the DOS and Z-factor results. 
There is a considerable reduction in velocity renormaliza¬ 
tion in the last three optical modes which indicates the 
polaronic feature of the carrier transports. 

FiglU shows the obtained total scattering rates 
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FIG. 5: (color online) Calculated scattering rates for total 
contribution(e-e and e-ph-e) and the e-ph-e contribution ver¬ 
sus energy. 


of iImJ2e-e+^'^'^e-ph-e) accordlug to the 
Matthiessen’s rule together with individual e-e in¬ 
termediated a phonon scattering rate 
where Im J2e-e '^e-ph-e denote scattering rates 

for Zr 02 at equilibrium structure and phonon-displaced 
snapshots, respectively. From the branches (d), (g), 
(h), and (i) of Fig 15] compared to the e-ph-e scattering 
rates in each vibronic modes, it is deduced that the 
contributions of the e-ph coupling play the main role in 
the total scattering rate in which cubic Zr 02 behaves 
metal-like having free carriers whereas, in panels (e) and 
(f) a weak contribution to the total scattering rates is 
illustrated, indicating a weak e-ph-e coupling in these 
two modes as we observed from the foregoing results. 
Moreover, there are small gaps and fluctuations in all 
configurations that can be attributed to the shape of 
the electronic band-structures. It is observed that in the 
higher conduction bands, there is a small band spread 
from point W to L while the band-width in the rest of 
the Brillouin zone, especially from F to X, experiences a 
wider energy variation, which is the cause of the decay 
in fluctuation rates. It is worth mentioning that our 
frozen phonon calculations are performed at the F point 
of the Brillouin zone where cubic Zr02 experiences its 
stable phase and does not show any soft modefsl fl5l|. 


Figl6| shows the obtained lifetimes for each phonon 
mode compared to the equilibrium state of Zr02- In 
acoustic branches, despite the electron assisted acous¬ 
tic phonon correlations, electron lifetimes remain unal¬ 
tered, which is in line with our foregoing results. In 
fact, in these modes electron lifetimes for inter-band and 
intra-band transitions have the same order of magnitude, 
though inter-band transitions due to the optical phonons 
are slightly larger than intra-band transitions due to the 
acoustic phonons. In optical branches, two effects are 
distinguished: 1) lifetimes in panels (d), (g), (h), and (i) 
decrease and follow the quadratic pattern of Liquid Fermi 
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FIG. 6: (color online) Electron lifetimes versus energy. Gon- 
trary to the optical modes, the calculated coupling between 
electron-acoustic phonons has no important effects on the 
electron lifetimes. The observed lifetime pattern in optical 
modes follow Fermi liquid model except for the (e) and (f) 
modes. 


model and 2) lifetimes in panels (e) and (f) increase. 
From the latter effect together our previous findings for 
the Z-factor and velocity renormalization, it is deduced 
that, the larger electron lifetimes result in the reduced 
scattering rates. Hence, the intra-band scattering is the 
more important transition process for these two modes. 
For the first effect, as it is stated from the Z-factor re¬ 
sults, the lifetimes of electrons are affected by the local¬ 
ized polaronic states consequently, in these optical modes 
inter-band scattering processes are relevant. Accordingly, 
to have low thermal conductivity Tph-ph (the phonon- 
phonon scattering life time which is an anharmonic ef¬ 
fect) should be small and Te_e (the bare e-e scattering 
life time) should be large. Referring to FiglU it is de¬ 
picted that a small Tph-ph is relevant with small Te-ph-e 
or equivalently, large values of which eluci¬ 

dated that the optical modes have the dominant part in 
low thermal conductivity of cubic zirconia. 


B. CONCLUSIONS 

In summary, we applied the GqWo@DFPT method 
within the frozen phonon model to analyze the phonon- 
assisted electron dynamics of cubic Zr02- We have 
traced phonons by observing the indirect electron- 
electron interaction to refrain from the time consum¬ 
ing calculations associate with the direct approach. De¬ 
spite the strong phonon-phonon interactions, from the 
obtained indirect electron-phonon coupling matrix along 
with the carrier velocity and energy renormalization, we 
observed a substantial polaronic character for all the ac¬ 
tivated lattice vibrations through the optical phonons ex¬ 
cept for the second TO and the first LO modes, making 
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these two dynamical configurations appropriate for prac¬ 
tical purposes. Furthermore, we observed a small to large 
polaron transition in the first LO mode. Our many-body 
GW calculations through the frozen phonon model, pave 
the way for the simple estimation of the order of electron- 
phonon and phonon-phonon interactions. 
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